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The stability of binary orbits can significantly shape the gravity wave signal which future Earth- 
based interferometers hope to detect. The inner most stable circular orbit has been of interest as 
it marks the transition from the late inspiral to final plunge. We consider purely relativistic orbits 
beyond the circular assumption. Homoclinic orbits are of particular importance to the question of 
stability as they lie on the boundary between dynamical stability and instability. We identify these, 
estimate their rate of energy loss to gravity waves, and compute their gravitational waveforms. 
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Coalescing compact objects are proficient emitters 
of gravity waves during the final minutes before their 
demise. They are the most popular candidate for the 
direct detection of gravitational radiation by the future 
LIGO/ VIRGO interferometers The gravitational 

wave frequencies sweep through the LIGO /VIRGO band- 
with of 1 — 10"* Hz as the inspiral gives way to the plunge 
and final coalescence. General interest in the onset of 
the instability to a plunge has focused attention on the 
innermost stable circular orbit (ISCO) The post- 
Newtonian (PN) approximation to general relativity ||] 
has been studied in detail and the ISCO determined in 
the absence of dissipation in Ref. 

There are other orbits which may be even more im- 
portant to the question of dynamical stability than the 
ISCO. The homoclinic orbit in particular arises at the 
boundary of the stable and unstable manifolds in the dy- 
namics. A homoclinic orbit begins and ends asymptoti- 
cally close to an unstable circular orbit, spiraling in and 
out around the center of mass. These trajectories define 
the region where dynamical instability, and even chaos, 
can develop [|j . Only black holes are compact enough to 
capture a companion in an orbit which approaches these 
smallest radii. Wc give a prescription for identifying ho- 
moclinic orbits in the 2PN expansion in the absence of 
dissipation, estimate the luminosity in gravity waves, and 
compute the gravitational waveform. 

Long lived binaries will likely circularize by the time 
of inspiral as angular momentum is lost to gravity waves. 
In dense stellar regions such as globular clusters the cap- 
ture of a companion can be an efficient mechanism for 
forming binaries There will not be sufficient time for 
these to circularize and we will have to accept them and 
their eccentricities. Recently some attention has been 
paid to gravity waves from more general orbits [Q. We 
use these advances to estimate the luminosity in grav- 
ity waves from a purely relativistic homoclinic orbit. In 
order to detect the waves, very precise theoretical tem- 
plates are needed. The raw data will then be compared 
to the template over many thousands of orbits to draw 
out a detection. To keep the data analysis computation- 
ally manageable, the theoretical templates are primarily 



for circular orbits. The homoclinic orbit will chirp up 
in frequency as it winds in toward the unstable circular 
orbit and then chirp down in frequency as it winds back 
out. To detect these, either more general templates, or 
new detection methods, are likely required. 



I. SCHWARZSCHILD BLACK HOLES 

Primarily for a point of comparison, we estimate the 
rate of energy lost to gravitational waves and the wave- 
forms for Schwarzschild black holes with a test mass com- 
panion in a relativistic orbit. In §|| we find the energy 
loss and waveforms in the 2PN expansion. Schwarzschild 
geodesic motion is confined to an orbital plane and is 
described by the familiar metric 
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There are two constants of the motion; the energy 
E and the angular momentum J per unit rest mass: 
(1 - 2^1) t' ^ E + I and K^cf)' = J, where a prime de- 
notes differentiation with respect to an affine parame- 
ter. In order to compare with the 2PN approximation 
we have used the binding energy E which is equal to 
the usual Schwarzschild energy minus one. Working in 
Schwarzschild time t, the conserved angular momentum 
is 
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and the energy constraint equation is 
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It is well known that for a large enough angular momen- 
tum, the black hole allows for two circular orbits; one 
unstable and one stable. The circular orbits are solu- 
tions to Vcs = 1. As was done in Ref. the location, 
energy, and angular momentum of the circular orbits can 
be parameterized as 
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As the angular momentum and energy drop the two cir- 
cular orbits move closer together and eventually meet at 
the ISCO. The ISCO is strictly speaking a saddle point 
of the dynamics and occurs for /3 = at R„„ = Rst = 
6, = V6/36 and J/m = 2^/3, E = 2\/2/3 - 1. 

The ISCO is of interest as it marks the onset of a dy- 
namical instability. There are other orbits which may 
be even more important to the question of the dynami- 
cal stability. The homoclinic orbits are defined as those 
trajectories which asymptotically approach the same un- 
stable orbit in the infinite past as they do in the infinite 
future. These orbits lie in the most unstable region of 
phase space and will be the site of any chaotic dynamics 
which may develop as a result of perturbations j^. In 
the binary system they begin infinitesimally close to an 
unstable circular orbit, wind out to a maximum radius 
and then wind back in again just reaching the unstable 
orbit. 

A trajectory which begins high on an unstable circular 
orbit has enough energy to escape the pair if its energy 
exceeds the asymptotic energy of a companion at rest 
at infinity; that is if i? > 0. The homoclinic orbit can 
therefore only exist if i? < at the unstable circular 
orbit. The first homoclinic orbit appears when /3 = 1/2 
and asymptotically reaches the unstable circular orbit at 
R*/to — A,m(j)^ = 1/8 and J^,/m = 4. Bombelli and 
Calzetta found the homoclinic orbit explicitly as 
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with the range < /3 < 1/2. 

No pair will ever reside exactly at the homoclinic or- 
bit. More naturally, we show a nearly homoclinic orbit 
obtained by starting a pair at the homoclinic J/m for 
/3 = 0.47 but with a radius 0.01m larger than the un- 
stable circular orbit. The trajectory is shown in fig. |^. 
The companion executes a few windings around the un- 
stable circular orbit before tracing the giant precessing 



ellipse. This path is much more interesting than a simple 
Keplerian elliptical orbit as is the gravitational waveform 
generated. 
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FIG. 1. = .47, R: 
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We calculate the luminosity in gravity waves and the 
waveform for the homoclinic orbits in the Newtonian 
quadrupole approximation. Strictly speaking, the New- 
tonian approximation is not appropriate for use with 
gravitational sources. In we utilize the waveforms 
properly estimated in the 2PN expansion of the two- 
body problem. The qualitative behaviour is strikingly 
similar to the quadrupole approximation of this section. 
The Schwarzschild solution may also have a role in the 
computation of the two-body expansion. In the test 
mass limit, the 2PN expansion converges slowly to the 
Schwarzschild solution. For this reason a hybrid expan- 
sion of the equations of motion has been proposed which 
is exactly Schwarzschild in the test mass limit Q . A sim- 
ilar hybrid expansion of the waveforms could reduce to 
these Schwarzschild expressions for a massive black hole 
with a light companion. 

The quadrupole moment is simply — mxiXj: 

I^^ — mR^ cos^ (j) 
lyy — mR^ sin^ 
Ixy — mR^ cos (j) sin (j) 

with all other components zero. The luminosity in grav- 
ity waves is estimated as 
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with the notation 



For a general orbit, this can be expressed as 
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with 

W ^^{r r +3rr) (1.9) 

Q = 2r R +6rr - 24RR02 „ i2r200 ^^q) 

P = 12rR(/) + 12r2(/) + 12RR0 + 2r2 -SR^f/)^ (1.11) 

The instantaneous rate of energy loss to gravity waves 
is shown in the top panel of fig. y. As expected, the lu- 
minosity is greatest when the binary is closest together 
although the pair spends a longer time at larger separa- 
tions. For 13 = 0, the homoclinic orbit and the unstable 
circular orbit are one and the same and the luminosity 
reduces to the expected circular orbit expr essio n. For 
the exact homoclinic orbit the solution eqn. 
substituted into eqn. 
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1.5) can be 



to write the luminosity as a 
function of (j) only. However, the expression is long and 
involved and not particularly illuminating. We present 
the numerical results instead. 

The transverse, traceless waveform can be estimated 
using the Newtonian quadrupole 
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with the Earth a distance D from the pair in a direction 
N, ^fj^ = Pijkrn^kni, and 

Pijkm = {Sik ~ NiNk) {6j„i - NjN,n) 

-liS,j-N,N,){6km-NkNra). (1.13) 

For the sake of illustration we place the Earth directly 
above the binary so that N = z. The waveforms for the 
two polarizations are = h^x = ^hyy 
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The observer can be moved without complication to an 
arbitrary location by using rotated quadrupole compo- 
nents, I = 10 with 0{N) the appropriate rotation 
matrix. A detector will respond to a linear superposition 
of /i+ and hx ■ The bottom panel of fig. |^ is the waveform 
h+ for the orbit of fig. 0. The middle panel focuses in on 
a shorter time interval to resolve the oscillations as the 
companion sweeps near the unstable circular orbit. This 
signal appears to chirp periodically. Although we have 
only used the Newtonian approximation, we will show in 
the next section that the 2PN estimation of the two body 
problem mimics these features with little alteration. 
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FIG. 2. Top: The luminosity in gravity waves for the orbit 
of fig. ^ as a function of t/m. Middle: The waveform hj^ over 
a short interval. Bottom: The waveform h+ over a longer 
interval. The scale of h+ is arbitrary. 

We show another nearly homoclinic orbit in fig. |^ for 
[3 = 0.3 which has the angular momentum of the unstable 
circular orbit but begins displaced by 0.01m from the 
circular radius. This orbit again winds tightly around 
the unstable circular orbit before tracing a wider orbit. 
The rate of energy loss to gravitational waves is shown in 
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the middle panel of fig. || and the waveform /i+ is shown 
in the bottom panel. 






FIG. 3. The orbit, waveform and the rate of energy loss for 
a nearly homoclinic orbit with j3 — .3,R = Run + 0.01, R = 0. 

This can be compared to the waveform for an elliptical 
orbit calculated in the Newtonian approximation [^,^. 
The elliptic orbits have similar peaks and drops but do 
not have the multiple oscillations inbetween. 



II. THE 2PN EXPANSION 

The Schwarzschild solutions of the previous section can 
be compared to the post-Newtonian approximations of 
the two-body problem. In the 2PN expansion, the center 
of mass equations of motion for the binary orbit can be 
written in harmonic coordinates as H,[KJ,[i1| 
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where m = mi + m2 is the total mass of the pair and for 
future use we define the product of reduced masses 77 = 
mim2/m^. In the test mass limit, 77 = 0, the harmonic 
coordinates are related to the familiar Schwarzschild co- 
ordinates by r = R — m. The form of A(r, r,^) and 
B{r^ r, (j)) depends on the order of the PN expansion and 
can be found in the Appendix. To 2PN order there are 
two constants of motion, the energy E(r, r, 0) and the 
angular momentum J(r,r,(j)) per unit rest mass. These 
elaborate expressions are also recorded in the Appendix. 
The two constants of motion ca n b e used in principle 
to eliminate r and (f> from eqn. (2.1). Just as with the 
Schwarzschild solution, the orbital dynamics can be re- 
duced to the study of one-dimensional motion in an ef- 
fective potential, r = —dVcft{r)/dr. Integrating we have 
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The constant is arbitrary and we shall hereafter set it to 
1. We will use eqn. (^) as a definition of Ves in the nu- 
merical solutions. A typical Ves for a black hole binary 
is shown in fig. ^ obtained by throwing the pair together 
at 2PN order and integrating to find Vcft = 1 — (l/2)f^. 
The shape of the potential does depend on both E and 
J/m but nonetheless this typical VeS gives an indication 
of the location of two circular orbits: one unstable circu- 
lar orbit near the maximum and one stable circular orbit 
near the minimum. 
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FIG. 4. The 2PN = 1 - (l/2)f2 with ri = 1/4 

for two trajectories. Top: Initial conditions are 
r/m = 20, r = -0.02, m<i> = 0.008. The potential depends 
on both the energy, E — —0.01298, and the angular mo- 
mentum, J/m = 3.3067. Bottom: The initial conditions 
are r/m = 6.8172, mcf) = 0.05176 which are ISCO values for 
77 = 1/4 but with r = 0.02 initially. 

To find the homoclinic orbit, we first find all of the 
circular orbits. For circular motion, f'o = fo = 0. From 



eqn. (2.1) this demands 
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We solve the circular orbit condition (2^) for (po{ro) as a 
function of the radius of the circular orbit and the ratio 
of masses rj. The solution is provided in the Appendix 
eqn. (A6) for reference and provides the complete initial 
conditions for circular orbits. We substitute (j)o into J to 
define an angular momentum along circular orbits 

• J„(r„) = J(r„,r„ = O,0„) (2.5) 

and into E to define an energy along circular orbits 

• Eoiro) = E{ro,ro = 0,^o)- (2.6) 

These are functions of only, once rj is specified. We 
find the smallest circular radius, r*, for which 



Eo{r,) < 0. 



(2.7) 



The first homoclinic orbit begins at this unstable circu- 
lar orbit. As the angular momentum and the energy de- 
crease, the last homoclinic orbit coincides with the ISCO. 
All homoclinic orbits can therefore be found by finding 
r» and rjgco s-s functions of rj. 
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FIG. 5. The energy along circular orbits (top) and the an- 
gular momentum along circular orbits (bottom) in the 2PN 
expansion with 77 — 0.1/(1.1)^. 

For a companion of mass 1/lOth its partner, rj = 
0.1/(1.1)^ then Eo{ro) and Jo(ro)/m are shown in fig. ||. 
For unstable circular orbits greater than r*/m = 4.2866, 
-E < 0. This can be repeated for any r\ to generate figs. 
^. These show the values of r^/m, the orbital frequency 
/* — (/)/(27r), and J^,/m for the unstable circular orbit 
when = as functions of rj. 
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FIG. 6. The radius (top), the orbital frequency (middle), 
and angular momentum (bottom) of the smallest unstable 
circular orbit for which there is a homoclinic orbit as functions 
of rj in the 2PN expansion. The energy by definition is zero. 

We follow Ref. g to locate the ISCO. The stability of 
the fixed points is tested by perturbing eqns. (2.1)-([2!^) 
about (ro, To = 0, (j)o) to obtain. 
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The perturbations obey 
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with iX the eigenvalues of (2.8), 

A = 0, A = ±(-a-6c)i/2. 



(2.9) 
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Stable oscillations about a circular orbit correspond to 
real A so that a + 6c < 0. The ISCO is the minimum 



value of To that has a real A and satisfies condition (2.4) 
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FIG. 7. along circular orbits in the 2PN expansion with 
r, = 1/4. 



Using eqn. (|Al|) for A and eqn. (A2) for B, we find 
(j)o{ro) in eqn. ( |A6| ). The resultant A^(ro) is illustrated 
in fig. (0) for an equal mass companion. For any 77, the 
values of risco, the orbital frequency /isco = 0/(27r), £^isco 
and Jisco are shown in figs. pM. 
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FIG. 8. The radius (top) and frequency (bottom) of the 
ISCO as functions of rj in the 2PN expansion. 
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FIG. 9. The energy (top) and the angular momentum (bot- 
tom) of the ISCO as functions of rj in the 2PN expansion. 

With the ISCO and the critical unstable circular orbit 
we have the initial conditions for any homoclinic orbit in 
the 2PN expansion. 

A. Nearly homoclinic paths 

An example of a nearly homoclinic orbit for rj = 
0.1/(1.1)2 obtained numerically by starting on an un- 
stable circular orbit at r/m — 4.409 and knocking the 
orbit slightly outwards. The binary orbit is not just a 
simple precessing ellipse. The orbit first glides around 
the unstable circular orbit, then executes a giant ellipse, 
swings around the inner unstable circular orbit again be- 
fore tracing another ellipse in a nearly opposite direction. 
This overall interleafed tracing then precesses. This orbit 
is relativistic and is not even gwasi-Keplerian. Snapshots 
of the orbit are shown in fig. |l^ and fig. |l^ shows the 
effective potential. A phase space diagram shows period- 
icity in (r, f) in fig. ^ 
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FIG. 10. A series of snap shots for a nearly homoclinic 
orbit. The ratio of masses is ni2/mi = 0.1 so that 
77 = 0.082644628. The initial conditions were r/m = 4.409, 
f = 0, = 0.0985406. Consequently E = -0.0084324135 
and J/m = 3.9166041778. 
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FIG. 11. The solid line is the lull potential while the dots 
indicate that part of the potential traced out by the relativis- 
tic orbit. 
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FIG. 13. The orbital frequency ol eqn. ( ^.12| ) measured in 
Hz as a function of t/m. 




FIG. 12. The phase space (r, r) for the trajectory of fig. |l^. 

The gravitational wave frequency measured by a dis- 
tant observer is roughly twice the orbital frequency for 
a circular orbit. For elliptical orbits the waves oscil- 
late at once, twice and three times the orbital frequency 
Since the orbital frequency here is not constant, the 
comparison is more ambiguous. The waveform oscillates 
roughly twice per winding around the center of mass. 
The orbital frequency in units of Hz is 



IAMq ) ^ ~ 2^ V L4GMq 



Hz, 



(2.12) 



and is shown in fig. |l^. For a pair with a total mass 
twice 1.4Mq, the frequency ~ 0(2000) Hz at the circular 
orbit and drops to 0(10) Hz. Although these frequencies 
are still within the LIGO range, the signal will appear 
to chirp up and then down again. For black hole pairs 
with larger total mass, the signal will move in and out of 
LIGO's bandwith. 



It is unclear if a circular orbit template integrated 
against data which chirps in and out will lock onto a 
signal. The far-zone energy and angular momentum flux 
is calculated in Ref. Q for general binary orbits in the 
2PN formalism. We use their results (quoted in the Ap- 
pendix) to draw Lqw in the top panel of fig. |l^. As 
expected the rate of energy loss is greatest near the un- 
stable circular orbit but the orbit spends much more time 
emitting gravity waves at the lower rate. 

Also shown is the waveform along this orbit using the 
2PN approximation to the far-zone waveform of Ref. . 

the 



We have writte n the results of Ref. 
Appendix eqn. (A13). Notice how strikingly similar the 
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2PN waveform is to the Newtonian approximation for the 
waveform from the similar Schwarzschild orbit of figs. |l| 
&0. 
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FIG. 14. Top: The rate of energy loss to gravitational 
waves for the orbit of fig. ^ as a function of t/m. Middle: 
Resolution of the waveform h+ for a short interval. Top: The 
waveform h+ over the entire interval of the orbit of fig. |l^. 
The scale is arbitrary. 

We find another nearly homoclinic orbit for equal mass 
binaries so that -q — 1/4. We begin the orbit at 



Ti/m = 5.01 m4>, = 0.0886. 



(2.13) 



These initial conditions are the same as the correspond- 
ing unstable circular orbit but with initial radius 0.01m 
greater. The orbit, luminosity in gravity waves, and the 
waveform h+ are shown in fig. |l^. Notice that the wave- 
form oscillates several times during the sweep around 
the unstable circular orbit. Again, there is a striking 
similarity between the 2PN waveforms and the wave- 
forms calculated in the quadrupole approximation for the 
Schwarzschild orbit of fig. |[ 





FIG. 15. Middle: Top: Anoth er nearly homoclinic orbit 
with initial values as in eqn. (2.13). The rate of energy loss to 
gravitational waves for the above orbit. Bottom: The wave- 
form h+. 

It is difficult to say how generic these relativistic orbits 
are in a realistic capture scenario. Clearly eccentric orbits 
are more likely than circular ones. For angular momenta 
in the range Jisco < < Jo, these relativistic orbits are 
most natural. If a black hole captures a companion in 
an initially eccentric quasi-Keplerian orbit. The angular 
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momentum and energy lost to gravity waves will not be 
sufficient to circularize prior to plunge but may drive the 
pair into the range E < and J < Jo and therefore onto 
a winding relativistic orbit. 

Although the frequencies are within the LIGO detec- 
tion bandwith, these waveforms may not be detectible 
by the method of matched filtering against circular tem- 
plates. Eccentric signals matched against circular filters 
result in a loss in signal-to-noise Q . The chirping up and 
down of these relativistic signal makes them even more 
difficult to pick up with a circular template. A quanti- 
tative analysis of the loss in signal-to-noise is currently 
underway. 

The effect of the radiative back reaction on the in- 
stability is also important. Dissipation due to radiative 
back reaction appears at 5/2PN-order. The location of 
the circular orbits is uneffec ted by these terms. The cir- 
cular orbits still obey eqn. ( |2.4D and A{ro, fo = 0, 0o) is 
the same with or without radiative back reaction since 
the 5/2PN contribution to A (eqn. ([A^)) is proportional 
to f. The stability of a given will however be affected. 
Dissipation will likely make orbits more unstable and in- 
tuition would suggest the ISCO would move outward. A 
binary on a nearly homoclinic path may plunge together 
as a result of energy loss which pushes the pair beyond 
the unstable asymptote. Alternatively, they may fall into 
a decaying oscillation about the stable point, gradually 
spiraling together. Which outcome results will depend 
on the initial conditions of the specific orbit. 

Even more intriguing is the possibility that the gravity 
waves disrupt the orbit to the point of chaos. A hypo- 
thetical gravitational wave can, under certain conditions, 
lead to a stochastic region along the homoclinic orbit . 
The importance of this is clear. Chaos may be an in- 
evitable element in both the coalescence and the gravity 
wave spectrum Q. However, the strength, frequency, 
and waveform specific to the binary may not have the 
right form to induce chaos. With the waveforms for gen- 
eral orbits of Ref. and the homoclinic orbits we have 
isolated, a more realistic search for chaos in coalescing 
binaries is possible. 
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APPENDIX A: 2PN EQUATIONS OF MOTION, LUMINOSITY AND WAVEFORMS 



The center of mass equations of motion for the relativistic binary can be written as 



r — ra 



{A + Br) 



TO „ r 



|JJlO|,|ll|. The equations for A, B, E and J can be found in Ref. |J,|10|,|ll[ and are included here for the purposes of 
keeping the paper self-contained. In the PN expansion, the expression for A and B can be obtained order by order in 
the small quantity v'^ — P + r'^cj)^. To 2PN order 



A = 1 + Ai + ^2 



(Al) 



with 
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-2(2 + + + 3?])v'^ - |r/f2 



-(12 + 297J) (^) + 7j{3 - ir^y + -rKl - 3r;)f 

^77(3 - 4ry)i;V2 - i?7(13 - Ari)—v^ 
2 2 r 

{2 + 25j] + 2ri^)—r^ 



and 



B — Bi + B2 



(A2) 



with 



Bi =-2(2-77)f 



1 . 
— r 
2 



77(15 + 4f])v'^ -{4 + 4I77 + 8f]'^) 377(3 + 277)7^2 



The subscripts refer to the order of the term in the PN expansion. To 2PN order, there is no dissipation yet included 
and there are two conserved quantities: The reduced energy and the reduced angular momentum. 
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3, ,1 1 , .o7n Ito, l/m\2 

+ - 1 - 377 t;^ + - 3 + 7^y- + -rj-r' + 

8 2 7'2r 2 \ r / 



E = -V 
2 

+ —(1 - 777 + lirf)v^ + -{21 - 2377 - 27r]^)—v^ + -tj(1 ~ 15r])—v^r^ 
16 8 r 4 r 

_ ^^(1 _ 377)-7^4 + i(14 - 5577 + 477^) f -)\^ + 1(4 + 6977 + 12,72) f -) \^ 
8 r8 Vr/8 \ r J 

. 3 



1 , N /Tn\ ' 

i(2 + 15.)(-) 



J = (r X v) |l iw2(l - 3,^) + (3 + + ^(1 _ 7,^ + 137/2)^4 

-I- 1(7 - 107/ - 977^)— - ^77(2 + bri)—r^ + 1(14 - 4I77 -I- 47;^) 
2 r 2 r 4 



TO 2 
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The energy and angular momentum are both measured in units of the reduced rest energy /i = 7712/771. 
For circular orbits ro = fo = Q which, from the equations of motion, requires 

• n mAn 



(A3) 



(A4) 



(A5) 



This condition reduces to a quadratic in 0^ and we obtain the angular frequency as a function of the radius of the 
circular orbit for 7; > 0: 
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'j>l[ro,v] 



F - VF^ - H 
K 



(A6) 



with 



F=\v{li-Av){-) +l-(l + 3,y)- 
2 \roJ To 

{ TTl \ ^ / Tn 3 /" TTt 

H = 477 (3 - 4?7) — 1 - 2(2 + 77)— + t(12 + 29??) — 



K^rl2ri (3-4??) 



m 



We show 4>t{To) as a function of the circular radius in fig. |T6|. This gives the complete initial conditions for any circular 
orbit, stable or unstable. 




FIG. 16. 4>l{ro) as a function of the radius of the circular orbit to 2PN order with 77 = 1/4. 
To include the effects of dissipation, the 5/2PN correction needs to be added to A and B. They are respectively 
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^-■q—r 
5 r 



5 r 



o 17 m 
3 r 



(A7) 
(A8) 



The instantaneous rate of energy loss to gravitational waves in the far-zone is calculated in Rcf. for general orbits 

/ \ inst 

to 2PN order. We cite their results here. Defining Lqw — ( 7^ I units of /i^, they derived 



far— zone 



LgW — Ejy + Eipjy + E2PN , 
8 r o T, 

i?iP^ = ^^{^[(785- 85277y 

- 2(1487- 1392r7)wV 

- 160(17-7;)— 



(A9) 
(AlO) 



+ 3(687 - 62077)?^^ + 8(36 7 - 157?)— 

r 



+ 16(1 - 477). 



(All) 
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2PN 
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2^„6 



42 



(1692 - 5497/7 + 4430?7^)i; 



— (17 1 9 - 102 7877 + 6292r]'^)v^r'^ 

-^(4446 - 523777 + I393772)— 
21 r 

^(2018 - 1520777 + 7572r]^)v^r* 
i(4987 - 8513?7 + 216577^)^ 

(281473 + 81828r; + 436877^) - 



(2501 - 2023477 + 8404?72)?^6 



756 
1 

42 

1 111 

— (33510 - 6O97I77 + 14290772)— ; 
63 r 

1 771^ 

— (106319 + 9798?7 + 5376772) — 
252 H 

2 

63 



^(-253 + 102677 - 5677^)^ 



(A12) 



The subscript denotes the order of the term in the PN expansion. 

The transverse traceless waveform to 2PN order was also computed in Ref. |^] and we write the instantaneous 
contribution here for completeness. The Earth is a distance D from the pair and in a direction N. The vector r is the 
vector marking the separation of the pair and n the corresponding unit vector. The mass difference is Sm — mi — 7712. 



The notation used is 



Tiirij and, as usual, n(^iVj^ — riiV 
2 



J -r 71^-7;, 

Sm 
m 



The waveform is then 



^..(1) + ^f..(l-5) , c^..(2) 



with the projection defined as 



(A13) 



and 



f..(0-5) _ 
SJJ — 



} 



( m 

{3(N.n)^ [277(,7;,) - 7^71,^] + (N.v) 

C^/'^ = ^ " ^'^^ [(N-n)'^ ((3«' - 15r2 + 7^) 77,, + 307^77(,7;,) - Uv, 
+ (N.n)(N.v)^ [l27^77,j - 3277(^7;^)] + (N.v)^ 

3(1 - 377)7;2 _ 2(2 - 37;)— J t-y + 4—7^(5 + 377)77(,7;j) 
3(1 - 377)7-2 _ (^iQ ^ 3^)^2 ^ 29— J I 

6/1-^) = —(1 - 277) ((N.n)3- [(457;2 - l05?-2 + 90-)7^77 
12 L r L r 
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427;2 _ 2107^2 ^ 8g_ ] ^^^^.^ 



(N.n)2(N.v) 



r 
, m 



21v^ - 1357^2 _^ 



- (N.n)(N.v)2— [487-77y - 18477(^7;^-)] 



77,>; + 3367^ni'i7; 
(N.v)^' 



1727;, 



777, 

4— 77y - 247;,j 



(A14) 
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^(N.n)^ I (69 - 66r?)t;^ - (15 - 90r])r'^ - (242 - 24r?)™ 

[(66 - 36r/)u^ + (138 + 84'n)r'^ 
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^(N.v) {[(23-1077) 
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+ (N.n)2(N.v)' 
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+ 
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(1242 - 417077 + 193077^)^ 
[1230 - 3I8O77 - 9077^] 7^7;^ } 

1 , ^ Q 777, 

60 r 
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(N.n)2— { [(117 - 48077 + 54077^)7;'* - (630 - 2850?/ + 4050772)7;27-2 
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777 
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(1284+ 1090??)- 



,21 



TO 



(915 - 1440?? + 720if)—f^ 
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[(132 + 54O77 - 580r;2)u2 + (300 - IMOry + 30077^)7^2 
(856 + 4007/ + 700772) — 
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TO 



The superscript denotes the order in the PN expansion. A fuU explanation of the derivation of these terms can be 
found in the original Reference . 
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